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SECTION  I 
INTRODUCTION 

V-GAMMA  Is  an  orbital  malfunction  simulation  computer  code  for  space  vehicles 
and  is  designed  for  use  in  assessing  the  range  of  atmospheric  reentry  conditions 
and/or  orbits  which  can  result  from  other  than  nominal  mission  trajectory  changes. 
It  gives  parametric  views  of  the  reentry  velocities  and  flight  path  angles  and 
orbital  lifetimes  which  may  occur  due  to  spacecraft  engine  operation  over  a  range 
of  conditions  of  initial  location  and  orientation,  and  engine  burn  duration. 

This  feature  extends  the  applications  of  the  code  to  Include  mission  design  studies 
and  optimization. 

The  calculational  methodologies  employed  are  presented  in  reference  1.  This 
technical  report  is  a  users  manual  and  no  attempt  is  made  to  represent  supporting 
equations  and  derivations.  However,  particular  sections  of  the  original  report 
are  referenced  frequently  as  a  source  of  more  in-depth  Information.  Appendix  A 
Is  a  cross  reference  of  variables  used  In  the  V-GAMMA  computer  code  to  the  vari¬ 
ables  employed  in  reference  1. 

Section  II  is  concerned  with  the  V-GAMMA  code  Itself  and  Section  III  is  con¬ 
cerned  with  SKETCH.  SKETCH  Is  a  display  code  written  to  digest  the  output  of 
V-GAMMA  into  a  series  of  maps  and  charts.  V-GAMMA  outputs  data  using  unformatted 
write  statements  onto  a  file  named  Tape  1  for  Input  to  SKETCH.  The  last  section 
presents  sample  problems.  Computer  output  for  the  sample  problems  is  given  in 
Appendix  B. 


1.  Balley-Souders,  A  Reentry  Environments  Velocity  and  Flight  Path  Angle  (V-GAMMA) 
Calculational  Methodology  for  (General  Space  Vehicle  Applications,  AfML-TR-74- 
244,  Air  Force  Weapons  Laboratory,  Kirtland  AfB  NM,  1976. 
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SECTION  II 
V-GAMMA  PROGRAM 

ACCIDENT  SCENARIOS 

Of  primary  concern  In  initiating  a  V-GAMMA  analysis  is  a  definition  of  the 
accident  scenarios  to  be  considered.  As  defined  in  Section  III  of  the  AFWL-TR- 
74-244  there  are  three  primary  effects  which,  in  general,  are  adequate  to  repre¬ 
sent  a  malfunction: 

Anomalous  time  of  burn  Early  burn 

No  burn 
Late  burn 
Nominal  burn 

Anomalous  burn  orientation  Random  orientation 

Random  within  a  specified  range 
Preferred  orientation 
Correct  orientation 

Anomalous  burn  duration  Nominal  burn 

Short  burn 

Extended  burn  or  to  depletion 

Any  of  the  above  effects  may  be  incorporated  individually  or  in  combination 
when  representing  a  malfunction.  Also  available  is  the  capability  to  treat  sub¬ 
sequent  burns  and/or  malfunctions.  After  a  specified  coast  period  any  one  or 
combination  of  the  above  effects  may  be  reincorporated  in  the  definition  of  a 
second  burn. 

When  simulating  accidents  with  specified  subsequent  burns,  the  code  will 
look  to  treat  secondary  burns  for  only  those  cases  which  result  in  trajectories 
with  perigees  greater  than  400  kft  altitude  and/or  drift  time  to  reentry  greater 
than  the  nominal  mission  coast  time  to  subsequent  burn.  Reentry  is  considered  to 
occur  at  400  kft. 

Accident  Probabilities 

After  the  accident  is  postulated,  one  can  then  formulate  its  probability 
of  occurrence.  In  general,  this  type  of  Information  Is  readily  available  for 
frequently  used  launch  systems.  However,  If  data  are  not  available,  some  type  of 
failure  modes  and  effects  study  may  be  warranted  Quite  often,  available 
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probabilities  are  presented  in  information  defining  idividual  system  component 
characteristics.  In  most  cases,  the  failure  of  any  one  of  several  components  may 
lead  to  a  particular  malfunction.  For  example,  the  failure  of  either  an  acceler¬ 
ometer  or  a  digital  computer  may  lead  to  a  randomly  oriented  thrust  vector. 
Therefore,  when  considering  a  malfunction  the  inverse  of  the  product  technique 
should  be  used  to  combine  the  individual  component  failure  probabilities.  On 
the  other  hand,  if  the  accident  being  postulated  is  a  combination  of  simulta¬ 
neous  malfunctions,  the  net  accident  probability  is  the  product  of  the  indivdual : 
malfunction  probabilities. 

The  code  takes  care  of  combining  accident  probabilities  for  sequential 
burn  simulations.  The  burns  are  considered  simultaneous  events  and  their  associ¬ 
ated  probabilities  are  appropriately  combined. 

Program  Logic 

The  calculational  methodology  and  logic  flow  employed  in  this  code  were 
given  in  Section  VIII  of  reference  1  and  will  not  be  reproduced.  However,  a  brief 
description  of  each  subroutine  and  the  main  program  is  provided.  Also  provided 
are  definitions  of  some  key  parameters  with  their  counterpart  symbols  as  employed 
in  reference  1. 

MAIN  PROGRAM— V-GAMMA 

The  bulk  of  the  program  logic  per  the  description  given  in  Section  VIII  of 
reference  1  is  located  in  the  main  program  (V-GAMMA).  The  program  is  well 
commented.  Loop  control  parameters  for  short  burns,  random  orientation  burns, 
and  multiple  burn  locations  can  be  readily  identified  and  related  to  the  V-GAMMA 
document. 

Input  and  output  control  statements  are  also  located  in  the  main  program. 
Accidents  are  inprocessed  in  block  form  by  a  call  for  subroutine  ACCIN.  ACCOUT 
prints  a  description  of  the  accident  to  be  treated  in  header  form.  This  appears 
before  the  list  of  possible  trajectories  which  could  result  from  the  burn  malfunc¬ 
tion.  In  general,  when  all  of  the  possible  cases  of  burn  time,  burn  location, 
orientation,  etc.,  are  considered,  a  multitude  of  possible  trajectories  will 
result.  Results  of  the  Individual  bums  are  stored  In  blocks  of  100  for  out- 
processing  as  a  string.  Those  trajectories  having  perigees  less  than  400  kft 
are  considered  reentries  and  their  descriptions  include  velocity,  angle  of  attack, 
location,  plus  other  critical  parameters  describing  the  reentry.  Those  trajec¬ 
tories  which  escape  reentry  will  be  either  orbital  decays  or  hyperbolic  or 
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parabolic  escapes.  In  this  situation,  a  description  of  the  resulting  trajectory 
is  given  with  the  reentry  parameters  set  equal  to  0.  Each  of  the  resulting  tra¬ 
jectories  is  assigned  a  probability  of  occurrence  derived  on  an  equal  basis  from 
the  overall  accident  probability  supplied  as  input  data. 

SUBROUTINES 

ACCIN 

ACCIN  reads  the  input  data  blocks  and  assigns  a  subscript  of  1  to  the 
parameters  describing  the  first  burn.  A  subscript  of  2  is  assigned  to  those 
parameters  describing  the  subsequent  burn  when  they  are  included  in  the  input 
data  list.  ACCIN  is  called  once  for  each  accident  and  the  entire  accident  is 
inprocessed  in  block  form. 

ACCOUT 

ACCOUT  is  called  immediately  after  ACCIN  to  print  the  data  just  inprocessed 
to  give  a  description  of  the  accident  being  treated. 

SELECT 

SELECT  Is  called  by  V-GAMMA  at  the  beginning  of  an  accident  simulation 
sequence.  Parameters  defining  the  accident  as  read  by  ACCIN  are  accessed  by 
SELECT  and  loaded  into  the  calculational  variables  used  by  the  main  program, 

V-GAMMA.  SELECT  is  recalled  if  the  accident  definition  has  indicated  a  subsequent 
burn.  In  this  situation,  subroutine  FREEZE  has  already  been  accessed  by  V-GAMMA 
to  store  information  necessary  to  restart  the  interation  through  the  possible 
parameter  variations  of  the  primary  burn  after  the  second  burn  is  specified  by 
the  Input  data  and  a  variation  of  the  primary  burn  does  not  result  in  a  reentry 
trajectory. 

UNIQUE 

For  accidents  of  random  oriented  thrust  vectors  within  some  predetermined 
solid  angle,  subroutine  UNIQUE  calculates  the  total  number  of  unique  orientations, 
NUMIS.  NUMIS  Is  used  to  calculate  equal  probabilities  for  the  various  orientations 
from  the  overall  accident  probability.  This  routine  Is  also  called  for  randomly 
oriented  second  burns. 

BURN 

BURN  performs  the  pulsed  bum  iteration  calculations.  The  Inputs  pro¬ 
vided  by  the  main  program  are  the  initial  orbit  and  the  vehicle  location  and  thrust 
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vector  orientation.  Also  provided  are  the  vehicle's  weight  before  and  after  the 
burn,  thrust  or  the  expected  change  in  velocity,  burn  time  and  the  number  of  burn 
pulses. 

Velocity  additions  are  simulated  using  pulsed  burns  rather  than  a  single 
change  in  velocity  at  the  point  of  burn  initiation.  By  approximating  the  burn  by 
a  number  of  pulses  with  appropriate  coast  periods  in  between,  the  actual  flight 
path  and  resulting  location  is  more  closely  approximated.  The  burn  simulation  is 
stopped  if  a  powered  reentry  occurs. 

RENTRY 

The  primary  function  of  RENTRY  is  to  set  the  values  of  the  variables 
in  the  output  list  describing  the  trajectory  resulting  from  a  burn  simulation. 

There  are  three  entry  points  to  the  subroutine.  They  are  BURNRE,  YESRE 

and  NORE. 

Entry  YESRE  is  accessed  when  it  has  been  determined  by  V-GAMMA  that  the 
perigee  of  a  calculated  trajectory  is  less  than  400  kft.  YESRE  computes  the  true 
anomaly  of  the  reentry,  its  velocity,  angle  of  attack,  and  other  critical  parameters. 

Entry  BURNRE  is  accessed  when  it  has  been  determined  that  reentry  occurred 
during  a  pulsed  burn  simulation.  It  computes  the  parameters  describing  the  reentry 
as  well  as  the  amount  of  fuel  remaining  on  board  and  the  remaining  bum  time. 

Entry  NORE  sets  the  output  parameters  describing  a  reentry  equal  to  0. 

This  condition  in  the  output  list  denotes  a  no  reentry  case. 

VECTOR 

VECTOR  is  a  multipurpose  vector  manipulation  routine.  It  has  three  entry 

points. 

Entry  VNORM— Given  vector  A,  VNORM  calculates  C,  the  magnitude  of  A, 
the  vector  B,  the  normal  of  A. 

Entry  V00T— -Given  vectors  A  and  B,  VDOT  calculates  C,  the  dot  product 

A'B. 

Entry  VCROSS— Given  vectors  A  and  B,  VCROSS  calculates  D,  the  cross 

product  A  X  B. 


7 


AFWL-TR-78-161 


NAXIS 

Given  the  radial  vector  R,  the  constant  C,  and  the  vector  K  the  unit 
vectors  I  and  J  which  form  the  orthonormal*  set  I,  J,  K  are- calculated  by  the 
relations 


T  =  R  -  C  K 

J  =  K  X  I 

NAXIS  is  employed  in  updating  the  vehicle's  heading  vector  with  changes  in  its 
location. 

VTRAN 

VTRAN  computes  the  elements  of  the  vector  X  in  the  equation 

X  (Transpose)*!!  *  Y  (Transpose)*T*B 
where  X  and  Y  are  vectors  of  length  3 

T  is  a  3x3  transformation*  natrix 
and  B  is  any  three  dimensional  orthonormal  basis 

STORE 

STORE  stores  relevant  information  describing  the  cases  resulting  from  an 
accident  simulation  for  output  as  a  string  by  subroutine  WDATA.  The  17  variables 
describing  a  case  are  stored  in  blocks  of  100. 

WDATA 

WDATA  outputs  information  placed  in  the  storage  arrays  by  subroutine  STORE. 
The  results  of  the  analysis  are  output  to  Tape  6  (output  file)  and  Tape  1  (data 
tape  used  as  input  to  program  SKETCH).  Data  is  placed  on  Tape  1  using  unformatted 
write  statements. 

STATUS 

STATUS  keeps  account  of  central  processor  time  consumption  and  provides 
this  information  at  the  end  of  the  output  list. 

LDRIFT 

LDRIFT  calculates  a  new  true  anomaly  Y  from  an  initial  position  X  after 
drifting  in  an  orbit  for  T  seconds. 
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TDRIFT 

Calculates  the  drift  time  in  seconds  required  to  drift  from  location  X 
to  location  Y  in  a  given  trajectory.  X  and  Y  are  in  radians.  The  drift  time,  T, 
is  always  positive.  The  true  anomalies  X  and  Y  determine  the  direction  of  travel. 
The  positions  supplied  to  this  routine  always  represent  clockwise  motion.  There¬ 
fore,  negative  drift  times  are  never  requested. 

F  (X ,E) 

Calculates  the  indefinite  integral  of  DX  /  (1 .+E*C0S(X) )**2  for  use  in 
trajectory  transit  time  calculations.  X  is  the  true  anomaly  in  radians  and  E  is 
eccentricity. 

FREEZE 

A  subroutine  for  storing  all  parameters  associated  with  a  variation  of 
the  main  burn  which  does  not  result  in  reentry  and  is  to  be  followed  by  a 
sequential  burn  simulation. 

Entry  START— This  part  of  the  subroutine  stores  all  parameters  defining 
the  first  burn  so  that  calculations  can  be  restarted  at  an  appropriate  location 
in  the  iteration  after  simulations  of  the  second  burn  have  been  completed.  It 
also  calculates  a  true  anomaly  for  the  second  burn  based  on  the  between  burn  coast 
period  supplied  as  input  data. 

Entry  STOP— This  portion  of  the  subroutine  is  called  after  all  iterations 
of  the  subsequent  simulation  have  been  completed.  It  restores  the  program  to  the 
conditions  which  existed  before  the  initiation  of  the  subsequent  burn  simulation. 

ZEROIN 

ZEROIN  is  called  by  LDRIFT  when  performing  trajectory  drift  calculations. 
ZEROIN  is  a  numerical  convergence  routine  and  is  used  in  Newton's  method  for  cal¬ 
culating  the  position  of  an  orbiting  or  escaping  vehicle  after  a  specified  coast 
period  from  some  initial  location.  Given  a  function,  F(X),  and  two  points,  AX  and 
BX,  such  that  F(AX)  and  F(BX)  have  opposite  signs,  ZEROIN  computes  an  X  between  AX 
and  BX  such  that  F(X)  =  0.0,  to  within  an  accuracy  specified  by  TOL. 

ACCIDENT  INPUT  DATA 

The  first  data  card  read  is  the  overall  title  of  the  analysis  being  performed. 

It  is  followed  by  a  card  giving  the  total  number  of  data  blocks  describing  accidents. 


9 


— ,mtt 


I 


AFWL-TR- 78-161 


The  data  blocks  follow  in  sequence.  Each  data  block  should  present  critical  burn 
parameters  completely  describing  the  bum  which  would  occur  as  a  result  of  the 
postulated  malfunctlon(s).  The  following  listing  describes  the  sequence  in  which 
the  data  cards  should  appear,  and  includes  the  card  number,  the  format,  and  the 
required  data.  Figure  1  shows  the  sequence  in  graphic  form  at  the  end  of  the 
listing. 


Card  No. 

Format 

Parameters 

1 

8A10 

Title 

Title  of  overall  analysis,  80 
alphanumeric  parameter  field 

2 

15 

NADB 

Total  number  of  accidents  in 
the  data  set. 

The  following  data  block  is  repeated  NADB  times. 

describing  each  accident  being 

treated: 

3 

215 

NACC 

Number  of  this  data  block 
in  the  NADB  block  sequence. 

NACC  identifies  the  data  blocks 
in  increasing  order  from  1  to 
NADB. 

NSUBB 

Identifies  this  accident  as 
having  a  subsequent  burn. 

This  version  of  V-GAMMA  can 
treat  only  one  subsequent  burn. 

So  NSUBB  will  be  0  or  1, 

4 

13A6.A2 

Title 

Title  of  the  individual  acci¬ 
dent,  80  alphanumeric  parameter 
field  . 

5 

F20.10 

PROBIN 

Probability  that  the  accident 
described  in  this  data  block 
will  occur, 

6 

2F10.5 

PERIN 

Perigee  of  the  initial  trajectory 

APOIN 

Apogee  of  the  initial  trajectory 
if  the  trajectory  is  an  ellipse. 
If  the  trajectory  is  other  than 
an  ellipse,  the  eccentricity 
should  be  given  instead.  (Alti¬ 
tudes  are  In  nautical  miles  above 
earth  radius)  • 

7 

3F10.5 

HVIN(3) 

Vehicle  heading  vector  at  loca¬ 
tion  of  burn.  The  heading 

vector  is  computed  from  the 
yaw  and  pitch  by  use  of  the 
following  transformation: 


10 


AFWL-TR-78-161 


8  2F10.5 


9  F10. 5 


10  215, 5X 

3F10.5 


11  2F10.5.15 


HVIN  -Hj  ■  sin  a  cos  b  0  0  I 

H2  *  0  sin  b  0  J 

Hj  s  0  0  cos  a  cos  b  K 

where  a  is  the  pitch  and  b  is 
the  yaw.  For  more  detailed 
description  see  section  V  of 
reference  1. 

W1 IN  Initial  weight  of  the  vehicle. 

W2IN  Final  weight  of  the  vehicle. 

The  difference  in  Ml  IN  and  W2IN 
represents  the  fuel  consumed  by 
the  engines  during  the  burn 
(pounds) . 

TCIN  If  a  sequential  burn  is  to  be 

treated  after  the  burn  being 
described,  TCIN  represents  the 
coast  time  between  burns. 
Otherwise,  TCIN  is  supplied  as 
0  (seconds). 

NBL  NBL  is  used  in  situations  where 

it  is  desirable  to  treat  a 
series  of  shorter  than  nominal 
burns.  If  NBL  is  specified  as 
greater  than  1 ,  then  the  burn 
specified  by  the  input  will  be 
treated  with  equal  probability 
of  having  NBL  different  burn 
lengths.  Specifically,  If  T 
is  indicated  bum  time  then  in 
the  Xth  iteration  the  burn  time 
is: 

Tx  =  T/NBL  *  X 
where 

1  <  X  <  NBL 

NIN  The  burn  specified  in  the  input 

is  divided  into  NIN  pulse  burn 
Iterations  to  more  closely 
aoDroximate  the  actual  path  of 
the  vehicle  during  the  burn. 

DTIN  The  net  burn  duration  is  seconds. 

THRUST  Vehicle  thrust  (lbs).  This 

parameter  is  not  necessary  if  a 
velocity  change,  DVIN,  has  been 
specified. 

TUN  True  anomaly  of  the  burn  (°). 

T2IN  For  situations  in  which  It  is 
desirable  to  treat  a  series  of 
possible  burn  locations.  T1IN 
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and  T2IN  specify  the  range  of 
true  anomalies  to  be  considered. 
Otherwise  set  T2IN  =  T1IN. 

NTIN  V-GAMMA  will  treat  NTIN  possible 

locations  within  the  T1IN  to 
T2IN  range  with  equal  probabil¬ 
ities  of  occurrence.  For  a 
single  burn  at  the  indicated 
location  set  NTIN  =  1. 

12  3F10.5  P1IN  &  Sets  the  upper  and  lower  limits, 

P2IN  respectively,  in  the  pitch  plane 

of  the  solid  angle  defining  the 
range  within  which  the  vehicles 
heading  vectors  can  vary  when 
simulating  a  series  of  random 
orientation  burns.  The  P1IN  and 
P2IN  angles  are  input  in  degrees 
relative  to  the  previously  speci¬ 
fied  heading  vector.  The  pitch 
plane  is  formed  by  the  radius 
vectors  and  vehicle  velocity 
vector. 

OPIN  The  random  orientations  will  be 

considered  to  occur  with  equal 
probability  in  increments  of 
DPIN  (°)  in  the  pitch  plane 
between  the  PUN  and  P2IN  angles. 

13  3F10.5  Y1 IN  Y1 IN  and  Y2IN  sets  the  maximum 

and  minimum  limits  in  the  yaw 
plane.  See  reference  1  for 
detailed  definitions  of  pitch 
and  yaw  angles. 

Y2IN  And  Y2IN  sets  the  limits  of  the 

solid  angle  defining  the  range 
within  which  the  vehicle's  head¬ 
ing  vectors  can  vary  when  simulat¬ 
ing  a  series  of  random  orientation 
burns.  Y1 IN  and  Y2IN  are  input  in 
degrees  relative  to  the  previously 
specified  heading  vector. 

DYIN  The  random  orientations  will  be 

considered  to  occur  with  equal 
probability  In  increments  of 
DYIN  (°)  in  the  yaw  plane 
between  the  Y1 IN  and  Y2IN  angles. 
Zero  yaw  specifies  the  pitch 
plane  of  the  heading  vector. 

The  following  cards  are  included  as  a  part  of  the  data  deck  if  NSUBB  is  equal  to 
1,  indicating  that  a  subsequent  burn  will  take  place. 
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14 

A80 

Title  (2) 

Title  of  subsequent  burn. 

15 

F20.10 

PR0BIN(2) 

Probability  that  burn  will 

occur. 

16 

2F10.5 

Ml  IN 

Weight  of  vehicle  before  burn. 

W2IN 

Weight  of  vehicle  after  burn. 

17 

215, 

NBL(2) 

See  data  card  10. 

5X, 

2F10.5 

IVIN 

See  data  card  10. 

OVIN 

See  data  card  10. 

DTIN 

See  data  card  10. 

18 

3F10.5 

P1IN 

See  data  card  12. 

P2IN 

See  data  card  12. 

DP  IN 

See  data  card  12. 

19 

3F10.5 

Y1  IN 

See  data  card  13. 

Y2IN 

See  data  card  13. 

DYIN 

See  data  card  13. 

20 

3F10.5 

HSEQ 

A  sequential  malfunction  can  assume  either  a  full  or  partial  failure  of  the  iner¬ 
tial  guidance  system.  This  read  statement  allows  the  flexibility  to  simulate  the 
sequential  malfunction  to  any  degree  of  inertial  guidance  failure.  If  HSEQ  is 
input  as  01  +  OJ  +  OK  subroutine  FREEZE  will  assume  a  complete  guidance  system 
failure  and  calculate  the  heading  vector,  assuming  a  spatially  fixed  heading  vec¬ 
tor  at  the  end  of  the  first  malfunction  burn.  If  HSEQ  is  inputted  other  than 
01  +  OJ  +  OK  this  inputted  value  will  be  used  as  the  heading  vector  in  the  sequen¬ 
tial  malfunction  calculations.  See  data  card  description  7  for  more  input  details. 
Note  that  if  random  orientations  are  specified  by  data  cards  15,  16,  and  17  the 
orientations  are  taken  about  the  subsequent  heading  vector  as  described  previously. 

The  above  described  data  cards  beginning  with  3  represent  one  data  block, 
here  should  be  NADS  data  blocks,  one  for  each  of  the  accidents  to  be  treated. 

V-GAMMA  OUTPUT 

As  discussed  previously,  many  malfunctions  can  result  in  a  multitude  of  possible 
trajectories.  The  V-GAMMA  output  consists  of  a  list  of  descriptions,  each  set 
describing  a  trajectory  which  could  result  from  the  specified  accident.  Each 
description  on  the  list  will  define  the  following  variales. 
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TITLE  OF  OVERALL  ANALYSIS 
NUMBER  OF  ACCIOENT  DATA  BLOCKS 


ACCIOENT  OATA  BLOCK 


SUBSEQUENT  BURN  DATA  BLOCK, 
IF  APPLICABLE 


ACCIDENT  DATA  BLOCK 


SUBSEQUENT  BURN  OATA  BLOCK, 
IF  APPLICABLE 


6/7 /8/9  CDC  DELIMETER 


Figure  1,  Sequence  of  Data  Cards 
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LEVEL 

THETA 

PITCH 


YAW 

P 

PR6 

E 

H 


T 

GEOP 

THETAC  - 
VC 

GAMMA 

ALPHAH  - 


BETAH 

TRE 

WFRE 

EVPROB  - 


Flag  Indicating  whether  description  refers  to  a 
primary  or  secondary  burn. 

Anomaly  at  which  malfunction  is  initiated. 

Initial  orientation  of  the  vehicle  in  the  pitch  plane 
formed  by  the  radius  vector  and  the  vehicle  velocity 
vector.  The  pitch  angle  is  relative  to  the  vehicle 
velocity  vector  ( degrees ) . 

Initial  vehicle  orientation  relative  to  the  pitch  plane. 

P  »  PRG*(l+E)Similatus  Rectum  of  the  trajectory  at 
burn  termination 

Perigee  of  the  resulting  trajectory 
Eccentricity  of  the  resulting  trajectory 

H  *  p  *  G  *  EM  Constant  of  angular  momentum  of  the 

trajectory 

EM  earth  mass 

G  universal  constant  of  gravitation 

True  anomaly  at  bum  termination 

Geocentric  perigee  of  the  resulting  trajectory 

True  anomaly  of  reentry 

Vehicle  velocity  at  point  of  reentry  (kft/s) 

Flight  path  angle  at  point  of  reentry  (degrees) 

Defines  vehicle  orientation  or  angle  of  attack  of 
reentry.  ALPHAH  is  the  angle  between  the  vehicle 
velocity  vector  and  its  heading  vector  and  lying  in 
the  orbital  plane  (degrees) . 

The  slip  angle  between  the  vehicle  heading  vector  and 
the  orbital  plane  (degrees). 

Time  (seconds)  from  malfunction  initiation  to  reentry 

The  amount  of  fuel  (lbs)  remaining  onboard  for  those 
cases  resulting  in  powered  reentry 

Probability  that  this  event  or  case  will  result  from 
the  defined  malfunction 
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SECTION  III 
SKETCH  DISPLAY  CODE 

SKETCH  Is  a  program  written  to  condense  the  output  of  V-GAMMA  into  a  set  of 
tables  and  charts  called  V-GAMMA  maps.  The  program  can  be  divided  into  two  parts: 
(1)  a  set  of  mapping  routines  which  generate  the  various  tables  and  charts  des¬ 
cribed  previously  and  (2)  a  set  of  routines  responsible  for  orbital  lifetime 
predictions.  These  two  parts  are  discussed  in  the  first  two  paragraphs  of  this 
section,  the  third  and  fourth  paragraphs  give  detailed  descriptions  of  input  and 
output.  The  capabilities  of  the  code  are  described  in  the  paragraph  defining 
the  input  parameters.  It  should  become  obvious  from  the  beginning  of  that  para¬ 
graph  that  SKETCH  will  allow  the  use**  considerable  amounts  of  flexibility  both 
in  selecting  accidents  for  review  from  the  overall  V-GAMMA  output  and  in  con¬ 
sidering  specific  cases  from  those  accidents. 

MAIN  PROGRAM  SKETCH 

SKETCH  contains  the  logic  necessary  to  make  decisions  and  call  appropriate 
subroutines  based  on  data  appearing  on  the  V-GAW1A  output  tape  and  the  specific 
V-GAMMA  maps,  charts,  and  tables  requested  by  the  input  data  cards. 

SUBROUTINES 

ZERO 

This  subroutine  zeros  the  various  storage  and  plotting  arrays.  It  Is 
called  initially  at  the  beginning  of  the  analysis  to  zero  all  a > rays  and  once 
after  each  accident  data  block  to  zero  those  which  are  used  to  generate  charts 
and  tables  that  are  specific  to  the  Individual  accidents.  Those  arrays  used  to 
generate  summary  Information  are  not  zeroed  after  the  initial  call  to  this 
subroutine. 

DECAY 

DECAY  is  called  by  SKETCH  when  an  orbital  trajectory  Is  being  processed. 
DECAY  In  turn  calls  those  subroutines  necessary  to  compute  an  orbital  lifetime 
for  the  trajectory.  As  discussed  in  the  first  paragraph  of  this  section,  two 
methods  are  available  for  orbital  decay  prediction.  For  long  lifetimes  (greater 
than  30  years)  a  three-dimensional  (3-0)  spline  interpolation  routine  is  used. 
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For  shorter  lifetimes,  L.  L.  Perini's  orbital  decay  program  is  used  directly 
(ref.  2).  DECAY  determines  which  method  should  be  used,  calls  the  appropriate 
subroutine,  then  calls  subroutine  PRINT  to  store  the  results. 

LIFE 

LIFE  is  an  orbital  lifetime  prediction  program.  The  program  performs  an 
orbit  by  orbit  iteration  decreasing  each  orbit  by  the  appropriate  amount  based 
on  atmospheric  drag.  In  keeping  track  of  the  elapsed  orbital  time,  LIFE  indexes 
KTO  (the  Zurich  support  number  indicator)  by  the  appropriate  amount  in  order  to 
take  into  account  changes  in  the  atmospheric  density. 

The  numerical  procedure  employed  in  the  analysis  is  based  on  the  develop¬ 
ment  presented  by  Klng-Hele  (ref.  3).  The  assumptions  employed  are  as  follows: 

1.  The  atmosphere  Is  spherically  symmetrical. 

2.  Air  density  varies  exponentially  with  height. 

3.  Lunar-solar  perturbations  are  neglected. 

4.  The  only  nonconservative  force  considered  is  air  drag  and  it 

acts  in  a  direction  tangent  to  the  orbital  plane;  i.e,,  normal  forces  are  neglected. 

5.  During  one  revolution,  the  action  of  air  drag  changes  the  orbit 
by  a  small  amount,  where  squares  can  be  neglected. 

6.  The  unperturbed  orbit  is  an  exact  ellipse. 

7.  The  atmosphere  rotates  with  constant  angular  velocity. 

8.  The  atmospheric  density  does  not  vary  with  time.  This  assump¬ 
tion  is  involved  in  the  derivation  of  the  equations  of  motion;  however,  long 
term  variations  are  considered  In  a  corrective  manner  by  reevaluation  of  the 
density  scale  height. 

ATM 

ATM  is  essentially  a  model  of  the  Jacchia  atmosphere  (ref,  4)  modified 
to  be  consistent  with  the  assumptions  given  in  the  description  of  subroutine  LIFE. 

2,  Perlnl,  L,  L,,  Orbital  Lifetime  Estimates,  John  Hopkins  Applied  Physics  Lab, 
ANSP-Mr8,  1974. 

3,  Klng-Hele,  D. ,  Theory  of  Satellite  Orbits  in  an  Atmosphere,  Butterworths ,  Inc. 
Wash  D  C,  1964. 

4,  Jacchia,  L.  G.,  "Revised  Static  Models  of  the  Thermosphere  and  Exosphere  with 
Empirical  Temperature  Profiles,"  Special  Report  SR-332,  Smithsonian  Astro- 
physical  Observatory,  1971. 
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Those  phenomena  which  can  cause  time  variation  in  atmospheric  density  are: 

1.  Variations  with  the  solar  cycle. 

2.  Variations  with  the  daily  change  in  activity  on  the  solar  disk. 

3.  The  diurnal  variation. 

4.  Variation  with  geomagnetic  activity. 

5.  The  semiannual  variation. 

6.  Seasonal-latitudinal  variations  of  the  lower  thermosphere. 

7.  Seasonal -latitudinal  variations  of  helium. 

8.  Rapid  density  fluctuations  probably  connected  with  gravity  waves. 

The  above  variations  are  in  general  short  term  with  the  exception  of 
those  due  to  changes  during  the  solar  cycle.  Short  term  effects  are  neglected 
in  this  model  because  of  computer  cime  requirements.  For  this  reason,  predictions 
shorter  than  10  years  are  not  ejected  to  be  very  accurate.  This  point  is  dis¬ 
cussed  in  reference  2  in  more  detail. 

ATM  also  calculates  and  fttfvss  modifications  to  the  density  scale  height 
(H).  This  parameter  is  used  to  account  for  changes  in  air  density  a  space  vehicle 
would  encounter  as  It  progress?*  through  an  elliptical  orbit. 

PRNTLF 

PRNTLF  is  called  by  subroutine  PRINT  once  for  every  100  orbital  lifetime 
predictions.  For  each  of  the  100  orbital  trajectories  subroutine  PRINT  prints 
the  Initial  vehicle  orientation  conditions,  the  critical  parameters  describing 
the  orbit,  the  predicted  orbital  lifetime,  and  the  probability  of  occurrence 
assigned  to  the  trajectory. 

STATUS 

This  subroutine  generates  a  summary  table  for  each  accident  data  block. 

The  tables  give  the  frequencies  and  probabilities  for  resulting  elliptical,  para¬ 
bolic,  and  hyperbolic  trajectories  with  breakouts  detailing  reentries,  escapes, 
and  orbital  decays.  If  THETAC  is  specified,  STATUS  summarizes  the  data  generated 
for  a  particular  theta  throughout  the  accident  block. 

M0DVG 

M0DVG  is  a  numerical  manipulation  routine  which  modifies  the  entries  of 
the  V-GAMMA  frequency  of  occurrence  maps  to  be  less  than  1000  for  plotting  under 
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13  format.  A  scale  factor  is  included  in  the  title  of  the  frequency  maps  which, 
when  applied  to  the  number  appearing  in  the  V-GAMMA  space,  yields  actual  frequency 
of  reentries. 

FMAP 

This  subroutine  generates  the  frequency  of  occurrence  V-GAMMA  maps. 

Similar  to  subroutine  PMAP,  subroutine  FMAP  plots  a  set  of  two  maps  for  each 
accident  simulation.  One  map  shows  all  reentries  and  the  other  shows  those  which 
are  powered.  A  set  of  overall  reentry  maps  displaying  all  reentries  resulting 
from  the  vairous  accidents  analyzed  is  provided  along  with  the  overall  probability 
maps. 

NUMEXP 

NUMEXP  is  a  numerical  manipulation  routine  which  evaluates  the  leading 
integer  (rounded  off)  and  power  of  10  exponent  from  reentry  probabilities  for 
entry  into  the  V-GAMMA  reentry  velocity  and  flight  path  angle  probability  maps. 

PMAP(K) 

This  subroutine  generates  the  probabilities  V-GAMMA  maps.  A  set  of  two 
maps  is  generated  for  each  accident.  One  map  shows  all  reentries.  Reentries 
which  occur  with  the  vehicle  engines  still  burning  are  also  displayed  in  a  second 
V-GAMMA  map  entitled:  Powered  Reentries,  this  second  map  will  not  appear  in  the 
output  list.  A  set  of  overall  V-GAMMA  maps  is  also  generated  displaying  all 
reentries  for  the  complete  set  of  accidents  analyzed. 

CHART 

CHART  generates  a  set  of  summary  tables  giving  the  results  of  each 
individual  accident  simulation.  It  is  called  by  SKETCH  as  data  are  processed 
from  the  V-GAMMA  output  tape.  There  are  several  entries  to  this  subroutine. 

ENTRIES 

ONE 

This  segment  of  the  subroutine  stores  data  into  the  matrix  for  the  fre¬ 
quency  charts  that  depict  the  time  to  reentry  from  malfunction  and  burn  time 
remaining  for  powered  reentry. 

TWO 

TWO  generates  reentry  time  charts  depicting  the  frequency  data  stored 
Into  the  reentry  time  matrix  by  ENTRY  ONE. 
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THREE 

This  segment  generates  the  matrix  that  is  used  to  produce  the  frequency 
charts  summarizing  the  vehicle  orientation  at  reentry  relative  to  its  velocity 
vector. 

FOUR 

This  entry  produces  the  charts  that  depict  the  frequency  of  reentry 
orientations  matrix  established  by  ENTRY  THREE. 

FIVE 

FIVE  sets  up  an  array  of  probabilities  versus  orbital  reentry  times. 

Data  are  simultaneously  stored  into  the  arrays.  One  contains  information  rela¬ 
tive  only  to  the  current  data  block  being  processed  from  the  V-GAMMA  output  tape 
by  SKETCH.  The  other  is  a  summary  storage  array  for  all  accidents  considered. 

SIX 

SIX  uses  the  array  established  by  ENTRY  FIVE  to  generate  an  orbital  life¬ 
time  probability  distribution  chart  for  each  accident  simulation. 

SEVEN 

SEVEN  is  called  by  SKETCH  just  prior  to  run  termination  to  generate  an 
overall  orbital  lifetime  distribution  chart  summarizing  the  results  of  all  of 
the  accidents  extracted  from  the  V-GAMMA  output  tape. 

SUBROUTINES 

PRINT 

PRINT  generates  frequency  plots  showing  the  distribution  of  orbital  life¬ 
times.  It  is  called  by  subroutine  DECAY  for  storage  of  each  orbital  lifetime 
calculated.  With  the  completion  of  each  accident  simulation,  SKETCH  calls  sub¬ 
routine  PRINT  for  the  resulting  distribution  of  orbital  lifetimes.  SKETCH  makes 
a  final  call  to  PRINT  just  prior  to  job  termination  for  a  combined  lifetime  dis¬ 
tribution  graph  of  all  accidents  simulated. 

SPIN I 3D 

SPLNI3D  is  called  by  subroutine  DECAY  for  lifetime  predictions  for  those 
orbits  expected  to  last  more  than  30  years.  SPLNI3D  interpolates  decay  time 
from  data  tables  located  near  the  beginning  of  SKETCH  as  follows: 
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Given  a  tabular  set  of  data 

(  (X(I) ,Y(J) ,2(1 fJ))  1=1 ... . ,NX  0*1 ... .  ,NY  ) 

that  describes  a  functional  relationship  between  the  independent  variables,  X  and 
Y,  and  the  dependent  variable,  Z,  SPLNI3D  uses  cubic  splines  to  interpolate  for 
the  value  of  the  function,  the  values  of  the  first  partial  derivatives,  and  the 
values  of  the  second  partial  derivatives  for  specified  values  of  the  independent 
variables. 

ORBITAL  LIFETIME  METHODOLOGY 

The  methodology  employed  in  SKETCH  for  predicting  orbital  lifetimes  is  based 
on  a  computer  model  developed  by  L.  L.  Perini. 

To  verify  the  accuracy  of  his  code,  Perini  ran  several  check  cases  for  a 
series  of  orbital  lifetimes  without  the  solar  effect,  using  tabular  values  of 

i 

the  1962  standard  atmosphere,  for  which  lifetime  solutions  are  available  based  on 
complete  numerical  solutions  (TMX-53385,  1966,  NASA).  The  resulting  comparison 
showed  excellent  agreement,  with  a  maximum  difference  of  about  8  percent. 

A  mean  solar  cycle  model  based  on  the  averages  of  cycles  8  to  19  was  incor¬ 
porated  in  the  code  (ref.  4).  The  results  led  to  the  following  conclusions: 

a.  With  the  mean  solar  cycle  as  a  basis,  lifetimes  are  three  to  four 
times  lonoer  than  would  be  oredicted  usinq  the  1962  standard  atmosDhere. 

b.  A  la  deviation  of  the  solar  cycle  from  the  mean  value  can  affect  the 
lifetime  of  a  factor  of  1.5. 

c.  For  lifetimes  greater  than  30  years,  estimates  based  on  a  nontime- 
varying,  reference  value  exospheric  temperature  derived  from  a  value  of  reference 
density  that  is  time-averaged  over  the  mean  solar  cycle  agree  well  with  more 
exact  calculations  that  directly  considered  the  mean  time-varying  solar  activity. 

d.  For  long  lifetimes  (>  30  yrs)  the  lifetime  is  directly  proportional 
to  the  ballistic  coefficient. 

e.  For  long  lifetimes  (>  30  yrs)  the  lifetime  estimate  is  independent 
of  where  on  the  solar  cycle  the  calculations  are  begun. 

Although  the  objective  of  this  work  was  to  provide  more  realistic  (rapid) 
estimates  of  orbital  lifetimes  by  accounting  for  solar  effects,  no  comparisons 
with  experimental  data  can  yet  be  made  because  of  the  long  lifetimes  under 
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consideration.  The  major  uncertainty  is  the  solar  cycle  model,  for  which  it  has 
been  assumed  implicitly  that  the  mean  of  future  cycles  will  be  represented  by 
the  mean  value  of  past  data. 

Based  on  the  above  generalizations,  the  lifetime  code  was  used  to  generate  a 
set  of  tables  of  Orbital  Lifetime/8allistic  Coefficient  versus  Perigee  and 
Eccentricity  for  use  in  a  cubic  spline  numerical  interpolation  and  differentia¬ 
tion  routine.  Predictions  of  orbital  lifetimes  for  the  extremely  high  altitude 
orbits  encountered  in  safety  analysis  were  too  time-consuming  for  efficient  use 
in  V-GAMMA.  However,  these  generalizations  are  only  applied  to  orbits  with 
expected  lifetimes  greater  than  30  years.  For  shorter  decay  times,  Perini's 
code  is  incorporated  directly. 

A  necessary  input  into  the  code  is  the  period  within  the  solar  cycle  at  which 
the  decay  is  to  be  initiated.  The  period  of  initiation  is  known  by  the  computer 
code  as  KTO.  A  graph  of  the  Zurich  sunspot  number  versus  time  into  the  solar 
cycle  is  included  for  the  purpose  of  evaluating  KTO  (figure  2).  Of  primary  con¬ 
cern  is  the  continuous  curve  which  represents  the  mean  of  cycles  8  to  19.  The 
dependence  of  the  lifetime  estimates  on  solar  activity  is  based  on  this  mean. 

Each  of  the  11  years  of  the  cycle  is  divided  into  quarters  so  that  KTO  can  range 
from  1  to  44  identifying  the  particular  quarter  of  the  year  in  which  the  calcu¬ 
lations  are  to  be  initiated.  KTO  is  evaluated  by  superimposing  the  current  11 
year  period  onto  the  graph.  This  is  done  simply  by  adding  some  multiple  of  11 
years  to  the  October  1964  date  such  that  the  date  that  the  decay  is  expected  to 
be  initiated  will  fall  into  the  11-year  period.  KTO  corresponding  to  the  quarter 
containing  that  date  is  obtained  directly. 

SKETCH  INPUT  DATA 

The  input  parameters  to  SKETCH  are  primarily  control  variables  for  the  various 
plotting  and  mapping  routines.  Through  the  correct  use  of  these  control  variables, 
the  data  tape  (Tape  1)  can  be  scanned  for  particular  data  sets  or  for  individual 
cases  of  interest.  The  exceptions  to  this  generalization  are  those  parameters 
employed  in  orbital  lifetime  predictions. 

Input  data  should  be  supplied  to  SKETCH  in  the  following  sequence  and  format: 
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Card  No. 

Format 

Parameters 

1 

(A80) 

Title  1 

The  overall  title  of  the  V-GAMMA  analysis 
80  alphanumeric  fields. 

2 

(15) 

KTO 

Index  for  Zurich  Sun  Spot  No.  See  sec¬ 
tion  on  Orbital  Lifetime  Methodology  for 
estimating  this  value. 

The  following  data  block  is  repeated  as  often  as  necessary  to  derive  the 
desired  summary  charts  and  maps  from  Tape  1. 

Card  No.  Format  Parameters 


1 

(A80) 

Title 

Title  of  individual  accident  or  case  to 
be  extracted  from  Tape  1 . 

2 

(F10. 4, 
215) 

BALL 

Ballistic  coefficient  of  satellite  con¬ 
figuration  after  the  malfunction. 

IC 

Sequential  accident  flag. 

IC  =  1 .  May  or  may  not  be  sequential 
malfunctions  on  the  data  files  of  Tape  1 
but  extract  only  the  results  of  primary 
bums. 

IC  3  2.  Definitely  sequential  burn  simu¬ 
lations  on  data  file  of  Tape  1  and  both 
will  be  extracted. 


3 


IB 


(15,  5X,  IFILE 
3F10.5) 


THETAC 


Level  of  accident  to  be  extracted  by  the 
present  data  set. 

IB  *  1.  Primary  accident. 

IB  =■  2.  Secondary  or  sequential  accident. 
NOTE:  Results  of  Primary  and  Secondary 
burn  simulations  for  a  single  accident 
should  be  extracted  in  sequence. 

The  number  of  the  file  on  Tape  1  which 
contains  the  data  set  to  be  treated  by 
this  data  block.  Each  accident  ir  entered 
on  Tape  1  as  a  separate  file.  T!v?  acci¬ 
dents  occur  in  the  same  sequence  that  the 
data  were  input  to  V-GAMMA.  The  output 
data  file  number  can  be  identified  as  the 
sequence  number  input  as  a  part  of  the 
V-GAMMA  input  data.  Data  resulting  from 
a  subsequent  burn  are  placed  on  the  same 
file  as  the  primary. 

The  maximum  value  of  true  anomaly  for 
selecting  cases  from  the  data  tape  (Tape 
1).  This  parameter  is  used  to  screen  for 
certain  anomalies  when  extracting  data 
from  Tape  1.  If  0  <  THETAC  <  360  only 
the  accidents  occurring  at  true  anomalies 
less  than  THETAC  are  extracted.  If  all 
cases  are  to  be  considered,  then  set 
THETAC  equal  to  370. 
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PITCHC  Only  cases  occurring  at  orientations  with 
pitch  equal  to  PITCHC  will  be  extracted 
from  the  data  tape.  If  all  pitch  orienta¬ 
tions  on  the  data  tape  are  to  be  extracted 
then  set  PITCHC  equal  to  370. 

YAWC  Only  cases  YAW  equal  to  YAWC  will  be 

extracted  from  the  data  tape.  If  all 
YAW  orientations  on  the  data  tape  are 
to  be  extracted  then  set  PITCHC  equal  to 
3f0. 


4 

(415) 

ITYPE(l) 

A  value  of  1  for  this  parameter  indicates 
to  SKETCH  that  frequency  V-GAMMA  maps 
should  be  generated. 

ITYPE(2) 

A  value  of  1  indicates  that  frequency 
V-GAMMA  maps  for  powered  reentries 
should  be  generated. 

ITYPE ( 3 ) 

A  value  of  1  indicates  that  probability 
V-GAf'WA  maps  should  be  generated. 

ITYPE(4) 

A  value  of  1  indicates  that  probability 
V-GAMMA  maps  should  be  generated  for 
powered  reentries. 

5 

FI  0.5 

BURNT 

Total  vehicle  burn  time.  This  parameter 
is  used  to  compute  the  remaining  burn 
time  for  powered  reentries. 

6 

6F10.5 

The  following  parameters  determine  the 

size  and  scale  range  of  the  V-GAMMA  maps. 
It  may  be  necessary  to  derive  the  scale 
limits  through  trial  and  error  in  order 
tr  obtain  maps  most  representative  of 
the  analysis. 

VMIN  -  Minimum  velocity  for  the  velocity 
scale  of  the  V-GAMMA  maps. 

VMAX  -  Maximum  velocity  for  the  velocity 
scale  of  the  V-GAMMA  maps. 

VCW  -  Channel  width  or  scale  increments 
for  the  velocity  scale. 

VMIN,  VMAX  and  VCW  should  be  selected  to 
give  a  maximum  of  24  channels.  This  num¬ 
ber  would  fill  a  standard  size  computer 
page  of  136  characters  horizontally.  The 
units  of  the  velocity  inputs  should  be 
kft/s. 

GMIN  -  Minimum  reentry  angle  for  the 

gamma  scale  of  the  V-GAMMA  maps. 

GMAX  -  Maximum  reentry  angle  for  the 

gamma  scale  of  the  V-GAMMA  maps. 

GCW  -  Channel  width  or  scale  increments 
for  the  gamma  scale. 
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GMIN,  GMAX  and  GCW  should  be  selected  for 
a  maximum  of  45  channels.  This  would 
vertically  fill  a  standard  size  printer 
page  of  64  lines  with  space  included  for 
a  title  and  scale  tables. 

SKETCH  OUTPUT 

The  V-GAMMA  maos  beoin  with  the  accident  definition.  The  accident  mav  be  due 
to  any  one  or  combination  of  the  following  malfunctions:  abnormal  length  of  burn; 
abnormal  time  of  burn,  and  misorientation  of  thrust  vector. 

Any  single  malfunction  or  combination  of  the  malfunctions  will  have  an  associ¬ 
ated  probability  of  occurrence.  It  is  often  desirable  to  specify  accidents  which 
may  occur  over  a  range  of  burn  lengths,  locations,  and/or  orientations.  In  this 
situation,  the  probability  of  occurrence  is  distributed  equally  over  the  specified 
range(s).  All  accident  iterations  which  result  in  a  reentry  at  a  given  velocity 
and  flight  path  angle  will  have  their  associated  probabilities  combined  and  placed 
in  the  indicated  V-GAMMA  block.  Probabilities  on  the  V-GAMMA  maps  are  read  as  a 
multiplier  times  an  exponent  of  10.  For  example,  2-6  is  read  as  2  x  10~6. 

Also,  generated  are  frequency  maps.  These  maps  indicate  the  total  number  of 
reentries  occurring  at  a  given  velocity  and  angle  of  attack.  Those  cases  resulting 
in  powered  reentries  are  extracted  from  the  overall  frequency  and  probability  maps 
and  displayed  in  a  separate  set.  In  addition  to  the  information  described  above, 
the  following  is  provided  as  a  set  of  summary  tables. 

a.  Maximum  and  minimum  velocity  and  anqle  of  attack  summary  of  points 
outside  of  the  range  selected  for  the  maps  (summary  of  off-scale  points). 

b.  Maximum  and  minimum  velocity  and  angle  of  attack  summary  of  on-scale 

points. 

c.  A  breakdown  of  the  quantity  and  types  of  trajectories  resulting  from 
the  accidents;  given  is  the  number  of  trajectories  which  fall  into  the  elliptic, 
parabolic,  and  hyperbolic  categories.  Those  trajectories  with  perigees  less  than 
400  kft  are  considered  reentries.  A  summary  is  given  for  each  trajectory  type 
telling  how  many  reentries  occurred,  and  of  those,  the  number  which  are  powered. 

It  should  be  noted  that  the  total  number  of  elliptical  escapes  actually  defines 
the  number  of  elliptical  orbits  with  perigees  greater  than  the  atmospheric  capture 
altitude  of  400  kft.  These  trajectories  result  in  orbital  decays. 
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d.  A  chart  giving  the  frequency  of  reentry  times  with  breakdowns  of 
time  from  malfunction  to  atmospheric  capture  and  remaining  burn  time  for  powered 
reentries. 

e.  A  chart  giving  the  frequency  of  vehicle  orientation  at  reentry  with 
a  summary  of  the  number  occurring  end-on,  side-on,  and  head-on. 

f.  A  chart  giving  the  probability  of  orbital  decay  times  with  breakdowns 
in  hours,  days,  months,  years,  and  centuries.  An  orbital  lifetime  distribution 
in  graphic  form  presents  probability  of  reentry  as  a  function  of  time  after  the 
accident.  This  latter  output  is  the  result  of  an  orbital  decay  model  developed 
by  L.  L.  Perini.  It  gives  good  orbital  lifetime  predictions  and  is  based  on  the 
1971  Jacchia  atmospheric  model. 

g.  A  set  of  overall  summary  maps  is  located  at  the  end  of  the  output 
list.  These  maps  represent  the  total  of  frequencies,  probabilities,  and  lifetime 
distributions  of  all  accidents  analyzed. 
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SECTION  IV 
SAMPLE  PROBLEMS 


V-GAMMA  SAMPLE  INPUT 

The  following  sample  problems  should  give  some  additional  insight  on  the  use 
and  applications  of  the  V-GAftIA  computer  code.  These  simulations  were  selected 
from  the  safety  analysis  performed  for  the  1977  Voyager- Jupiter-Saturn  fly-by 
mission.  The  nominal  mission  profile  called  for  two  engine  burns  after  an 
initial  orbit  was  achieved.  The  first  was  a  second  burn  of  the  Centaur  liquid 
propellant  booster,  and  the  second  was  a  burn  to  depletion  of  the  solid  propel¬ 
lant  propulsion  module. 

The  accident  scenarios  present  a  description  of  the  accident  along  with  the 
vehicle  burn  and  trajectory  parameters,  followed  by  listing  (figure  3)  showing 
the  data  cards  required  to  input  the  described  accident.  The  V-GAMMA  output 
for  the  described  accident  is  presented  in  appendix  B.  A  detailed  description 
of  the  output  which  the  code  generates  is  given  in  the  paragraph  entitled 
V-GAMMA  Output  (section  II). 

The  output  generated  by  V-GAMMA  and  SKETCH  for  the  sample  problems  is  given 
in  appendix  B.  A  system  dayfile  listing  is  also  in  appendix  B.  It  is  included 
as  an  aid  to  the  user  in  setting  up  the  V-GAMMA  and  SKETCH  execution  sequence. 

SHORT  BURN 

This  accident  presents  a  situation  in  which  there  is  a  premature  engine  shut¬ 
down.  The  key  information  needed  for  this  simulation  is  the  nominal  engine  burn 
time  and  the  number  of  burn  durations  initiating  at  the  nominal  burn  location  to 
be  considered.  The  various  burn  durations,  which  range  within  the  nominal  burn 
time,  are  treated  with  equal  probability  of  occurrence.  The  input  parameters 
include  a  nominal  burn  time  of  356.51  seconds  with  an  initial  and  final  engine 
weight  of  33150.362  and  9404.316  pounds,  respectively.  Engine  thrust  is  29555.0 
pounds.  Engine  thrust  is  simulated  using  15  pulse  burn  iterations.  The  malfunc¬ 
tion  simulation  treats  five  burn  lengths  with  a  conditional  occurrence  probability 
of  0.3333.  The  burns  are  initiated  at  a  true  anomaly  of  8.32321°  in  an  elliptical 
orbit  having  a  perigee  of  83.363  nmi  and  an  apogee  of  101.92  nmi  with  a  vehicle 
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V0YA6CR  SAMPLE  V-GAMMA  PROBLEM 
3 

1  0 

SHORT  ON  Time  correct  ORIENTATION  BURN  (SECOND  centaur  »urn> 

.  3333 

93.363  101.92 

-.1442  .02749  .98917 

331*0.362  9404.3160 
1*14.509 

5  25  17422.313  356.51 

3.23321  8.23321  1 

0 .  0 .  0 . 

0 .  0.  0 . 

2  1 

NOMINAL  ON  TIME  RANDOM  ORIENTATION  BURN  (SECOND  CENTAUR  BURN) 

.3333 

83.363  101.92 

0.0  0*0  1.0 

331*0.362  9404.3160 
184.509 

1  25  356151  29555. 

8i23321  8.23321  1 

0.  360.  15. 

0.0  90.  45. 

SUBSEQUENT  ON  TIME  MISORIENTEO  PROPULSION  MODULE  BURN  (RAN80M  CENTAUR) 
.999976 

4438.34  2100.67 

1  OS  44*7992  14633.229 

0.  12.  12. 

0.  25.  25. 

0.0  0.0  0.0 

3  1 

NOMINAL  BURN  AT  A8N0RMAL  LOCATION  (SECOND  CENTAUR  BURN) 

;33J3 

83.363  101.92 

-.1442  .02749  .98917 

331*0:362  9404.3160 
184*509 

1  25  356:51  29555. 

8.23321  33.23321  5 

0.0  0.0  0.0 

0.0  0.0  0.0 

RANDOMLY  ORIENTED  PROPULSION  MOOULE  BURN  (ABNORMAL  CENTAUR  BURN  LOCATION! 
1080024 

4438:34  2100.67 

1  05  6340.89  44.7992 

-43.  45.  15. 

0 .  0 .  0 . 

0.0  040  0.0 

0000000000009000000000 

Figure  3.  V-GAMMA  Sample  Input 
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heading  vector  of  -  0.1442  I  +  0.02749J  +  0.98917K.  The  burn  periods  are  59.42, 
118.83,  178.25,  237.67  and  297.09  s. 

With  the  exception  of  the  burn  duration,  all  other  mission  parameters  are 
nominal . 

RANDOM  ORIENTATION  BURN  WITH  MISORIENTED  SUBSEQUENT  BURN 

The  vehicle  and  initial  conditions  described  for  the  short  burn  simulation  are 
employed  for  the  first  burn  of  this  simulation.  However,  inst*ead  of  a  premature 
engine  shutdown,  a  guidance  system  failure  is  treated.  The  random  orientation  is 
simulated  by  specifying  pitch  and  yaw  variations  of  O’ to  360°  and  0.0°  to  +  90°, 
respectively,  about  the  vehicle’s  heading  vectors  to  define  a  hemispherical  solid 
angle.  This  takes  advantage  of  symmetry  of  events  with  respect  to  plus  and  minus 
yaw  angles  to  represent  complete  spherical  variations  in  orientation.  Pitch  and 
yaw  increments  of  45°  specify  the  angular  distances  between  the  vehicle's  heading 
vectors  at  which  the  various  burns  are  initiated.  A  conditional  occurrence  proba¬ 
bility  of  0.3333  is  divided  equally  among  the  equally  spaced  heading  vectors  by 
V-GAMMA.  For  each  of  the  heading  vectors,  the  vehicle  burn  is  simulated  followed 
by  a  coast  period  of  184.509  s.  After  the  coast  period,  a  second  burn  is  initiated. 

In  setting  up  a  burn  simulation  in  which  a  second  burn  follows  a  burn  during 
which  a  malfunction  was  experienced,  an  assessment  should  be  made  as  to  how  the 
malfunction  may  effect  the  second  burn.  Two  scenarios  which  may  be  descriptive 
of  the  second  burn  conditions  are: 

a.  A  complete  guidance  system  failure  resulting  in  the  vehicle's  main¬ 
taining  the  inertial  heading  of  the  first  burn  during  the  coast  period  and  subse¬ 
quent  second  burn. 

b.  The  vehicle  attempts  to  assume  the  nominal  heading  for  the  second 
burn  under  conditions  in  which  it  is  unaware  of  its  heading  or  location  due  to 
the  first  malfunction. 

For  instructive  purposes  the  second  scenario  will  be  treated  by  assuming  two 
burn  orientations  both  having  equal  probability  of  occurrence.  The  first  head¬ 
ing  will  be  the  inertial  heading  of  the  vehicle  at  the  time  of  the  first  burn 
and  the  second  will  represent  a  specified  pitch  and  yaw  rotation  from  that  heading. 

The  second  burn  has  an  associated  probability  of  0.99997.  It  will  be  simulated 
using  five  pulse  burn  iterations.  The  final  stage  has  a  net  thrust  of  14633.229 
pounds  and  has  initial  and  final  weights  of  4438.34  and  2100.67  pounds,  respectively. 
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The  simulation  is  set  up  to  treat  two  burns,  one  along  the  inertial  heading 
which  the  vehicle  has  at  time  of  first  burn  and  a  second  at  a  yaw  of  25°  and  a 
pitch  of  12°  relative  to  that  heading.  V-GAMMA  divides  the  above  probability 
of  occurrence  by  2  so  that  the  two  conditions  of  the  second  burn  are  considered 
to  have  equal  probabilities  of  occurrence  of  0.499985. 

ABNORMAL  BURN  LOCATION  WITH  SUBSEQUENT  MISDIRECTED  SECOND  BURN 

This  simulation  is  designed  to  represent  a  situation  in  which  the  vehicles 
burn  takes  place  at  an  abnormal  true  anomaly.  The  first  burn  is  followed  by  a 
propulsion  module  burn  which  has  a  random  pitch  orientation.  As  with  the  random 
orientation  burn,  this  simulation  is  performed  by  specifying  a  range  over  which 
the  burn  specifications  can  vary.  With  the  exception  of  the  location  of  the 
first  burn  and  the  pitch  orientation  of  the  second  burn,  all  mission  parameters 
are  normal.  The  abnormal  burn  location  probability  is  0.3333.  It  is  distributed 
over  a  range  of  45°  beginning  at  the  nominal  true  anomaly  of  8.323°.  A  frequency 
of  five  burn  locations  within  this  range  is  treated.  After  separation  of  the 
coast,  the  propulsion  module  burn  is  simulated  with  a  random  pitch  orientation 
about  the  nominal  mission  heading  vector.  The  heading  vector  for  the  second  bum 
is  specified  in  the  input  as  0.314311  2.08  x  10“ 5  J  and  0.949322  K.  This  vector 
was  derived  from  pitch  and  yaw  angles  of  18.319  and  0.0012,  respectively,  using 
the  transformation  matrix: 


H»“ 

= 

sin  a  cos  b  0 

0 

= 

0  sin  b 

0 

j 

»3 

a 

0  0 

cos  a  cos  b 

K 

- 

- 

- 

where  a  is  the  pitch  angle  and  b  is  yaw.  The  vehicle  heading  at  second  burn 
initiation  Is  varied  within  the  pitch  plane  about  the  nominal  heading.  Seven 
burns  are  simulated  over  a  pitch  range  of  +45°,  -  45°,  about  the  nominal  head¬ 
ing  Increments  of  15°.  The  probability  that  the  second  burn  will  occur  at  an 
abnormal  pitch  orientation  is  0.00002.  The  option  of  specifying  velocity  change 
Instead  of  engine  thrust  is  employed  for  this  burn.  The  nominal  velocity  change 
is  6340.89  ft/s. 
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SKETCH  SAMPLE  INPUT 

The  SKETCH  sample  problems  are  intended  to  illustrate  the  flexibility 
offered  in  choosing  data  from  V-GAMMA  output  tapes  to  be  presented  in  chart-form 
and  in  identifying  the  types  of  charts  to  be  generated.  The  problem  analyzes 
the  data  output  tape  of  the  previously  discussed  V-GAMMA  sample  simulations. 

For  instructive  purposes,  only  portions  of  the  data  resulting  from  the  various 
accident  cases  are  extracted  from  the  tape.  In  addition,  the  accidents  are 
extracted  in  a  sequence  different  from  the  sequence  in  which  the  accidents  were 
simulated  using  V-GAMMA  and  placed  on  the  tape.  Only  four  burn  simulations  are 
considered.  Since  a  data  set  is  required  for  each  reading  of  a  file  from  the 
V-GAMMA  output  tape,  four  data  sets  follow  the  first  two  input  cards.  If  a 
single  file  is  to  be  read  more  than  once,  then  an  appropriate  number  of  sets  are 
required  for  that  file.  Figure  4  is  a  listing  of  the  SKETCH  sample  input  data 
cards . 

Title  of  the  Simulation 

The  first  data  card  of  the  SKETCH  input  data  gives  an  overall  title  for 
the  accident  simulations.  The  second  card  gives  a  value  of  9  for  the  KTO  input 
parameter.  This  indicates  that  a  starting  point  in  the  11  year  solar  cycle 
appropriate  for  January  1978  should  be  used  for  orbital  lifetime  calculations. 

Subsequent  Misoriented  Propulsion  Module  Bums  following  the  Random 
Orientation  Centaur  Burns 

All  information  giving  the  results  of  the  subsequent  misoriented  propul¬ 
sion  module  burns  which  followed  the  random  orientation  centaur  burn  simulations 
are  presented  in  V-GAMMA  maps.  The  centaur  burn  simulations  which  resulted  in 
reentries  were  not  followed  by  propulsion  module  burns  and,  therefore,  do  not 
appear  in  the  maps.  This  first  data  card  for  this  set  gives  the  title  of  the 
data  file  to  be  analyzed,  A  ballistic  coefficient  of  0.2106  units  is  used  for 
the  propulsion  module.  It,  along  with  a  sequential  burn  flag  of  2  Indicating 
that  the  data  file  contains  sequential  burns  and  a  level  number  of  2  indicating 
that  the  subsequent  burns  should  be  extracted  appear  on  the  second  data  card. 

The  third  data  card  gives  a  file  number  of  2  for  the  data  set.  This  SKETCH 
tells  that  the  second  file  of  the  V-GAMMA  tape  contains  the  information  giving 
the  results  of  the  simulation  being  considered.  In  addition,  this  card  gives 
THETAC  PITCHC  and  YAWC  values  of  370°  indicating  that  SKETCH  should  extract  all 
cases  of  true  anomaly,  pitch  and  yaw.  The  fourth  data  card  gives  ITYPE  (1), 
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ITYPE  (2),  ITYPE  (3),  and  I TYPE  (4)  all  values  of  1.  This  indicates  that  fre¬ 
quency  and  probability  V-GAMMA  maps  should  be  generated  for  both  free  flight 
and  powered  reentries.  An  engine  burn  time  of  44.7992  s  is  placed  on  the  fifth 
data  card.  The  final  data  card  for  the  set  identifies  the  scales  to  be  used  for 
the  V-GAMMA  maps.  Indicated  from  left  to  right  are  minimum  and  maximum  veloci¬ 
ties  of  14.0  and  38.0  ft/s  with  a  scale  increment  1.0  kf/s  and  minimum  and  maxi¬ 
mum  flight  path  angles  of  0.0°  and  90.0°  with  a  scale  increment  of  2.0°. 

Abnormal  Location  of  Centaur  Burn  and  Subsequent  Propulsion  Module  Burn 

Both  the  first  and  second  burns  of  this  accident  simulation  are  presented 
in  the  V-GAMMA  maps.  Individual  sets  of  maps  are  generated.  The  first  set  gives 
the  results  of  the  Centaur  burn.  Simulations  and  the  second  set  gives  those  of 
the  subsequent  propulsion  module  burns.  A  set  of  SKETCH  input  data  cards,  simi¬ 
lar  to  that  for  the  subsequent  misoriented  propulsion  module  burn,  is  required 
for  each  of  the  two  burns.  The  data  card  following  the  title  card  for  the  abnor¬ 
mal  burn  location  gives  a  ballistic  coefficient  for  the  Centiur  vehicle  of  0.00395. 
Also  on  this  card  is  a  sequential  burn  flag  of  2  and  a  level  number  of  1  indicating 
that  simulations  of  the  first  burn  are  to  be  extracted  for  this  data  set.  With 
the  exception  of  a  file  number  of  3  and  a  burn  time  of  363.1347  s  for  the  Centaur, 
the  remainder  of  the  data  set  is  identical  to  that  of  the  misoriented  propulsion 
module  burn.  The  data  set  for  the  propulsion  module  burns  following  the  abnormal 
location  burns  of  the  Centaur  is  identical  to  the  data  set  for  the  misoriented 
propulsion  module  burn  with  the  exceptions  of  a  file  number  of  3  and  a  value  of  0 
for  both  ITYPE  (1)  and  ITYPE  (2)  indicating  that  frequency  and  probability  maps 
are  not  to  be  generated  for  powered  reentries. 

Short  on  Time  Correct  Orientation  Centaur  Burn 

Data  resulting  from  this  simulation  are  located  on  the  first  file  of  the 
V-GAMMA  output  data  tape.  The  sequential  burn  flag  is  set  to  1  with  a  level 
number  of  1  indicating  that  no  subsequent  burn  simulations  are  to  be  extracted 
from  the  data  file  and  that  orbital  lifetimes  should  be  calculated  for  the 
results  of  the  first  burn.  The  file  number  in  the  SKETCH  input  data  set  is 
appropriately  given  a  value  of  1.  The  remainder  of  the  data  set  is  identical 
to  that  for  the  simulation  of  the  abnormal  Centaur  burn  locations. 


34 


4 


AFWL-TR-78-161 


APPENDIX  A 

CROSS  REFERENCE  OF  VARIABLES 

FOR  THE  V- GAMMA  COMPUTER  CODE 
ANO  THE  V-GAMMA  METHODOLOGY  REPORT  (ref.  I) 
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V-GAMMA  REPORT* 

COMPUTER  PROGRAM 

“  (0) 

alpref 

- 

Definition  of  reference 
trajectory 

■ 

ALPHA 

- 

Interim  value  during  pulse 
burn  calculations 

aH 

ALPHAH 

- 

Calculated  in  RENTRY  60 

SH  *(pp  42) 

BETAH 

- 

Calculated  in  RENTRY  63 

C  (pp  24, Eq  20) 

C 

- 

Calculated  in  VECTOR  21 

Y 

Gamma 

- 

Calculated  in  RENTRY  38 

e 

EREF 

Definition  of  reference 
trajectory.  Calculated  in 
SELECT  33 

E 

- 

Initial  interim  value  dur¬ 
ing  pulse  burn  calculations 

EP 

Interim  value  in  primed 
coordinate  system  during 
pulse  burn  calculation. 
Calculated  in  BURN  81 

h 

HREF 

* 

Definition  of  refernce  tra¬ 
jectory  calculated  in 

SELECT  39 

H 

- 

Initial  interim  value  during 
pulsed  burn  calculation 

HP 

Interim  value  in  primed 
coordinate  system  during 
pulse  burn.  Calculated  in 
BURN  77. 

H  (0) 

HVIN(I) 

- 

Input  data  ACCIN  40 

L 

NTIN 

- 

Input  data  ACCIN  51 

m  (pp  32, Eq  47) 

EOVERM 

- 

Calculated  in  BURN  80 

m 

Used  to 

calculate  g  for  depletion  burns 

This  is 

not 

done  in  V-GAMMA 

M 

TUN 

- 

Input  parameters  ACCIN  51 

THETA 

- 

Initial  vehicle  location 
during  V-GAMMA  calculations 

T 

- 

Initial  vehicle  location 
during  burn  iteration 

♦Reference  1,  pages  66-69  and  cited  pages. 


AFWL-TR-78-161 


n,  N 

P 


P 


TIP 

T2P 

NIN(I) 

N 

PREF 

P 

PP 


Initial  vehicle  location 
after  burn  phase  of  burn 
iteration 

Final  vehicle  location  after 
coast  phase  of  burn  iteration 

Input  parameter  ACCIN  47 

Do  loop  parameter  for  burn 
iteration  BURN  49 

Definition  of  reference 
trajectory  Calculate  in 
SELECT  39 

Initial  interim  value  during 
pulse  burn  calculations 

Interim  value  in  primed 
coordinate  system  during 
pulse  burn  Calculated  in 
BURN  79 


These  angles  are  determined  from  the 
inputted  heading  vector  and  are  used  to 
establish  the 


n 

n  n  coordinate  system 

e 

ez  e» 

V-GAMMA 

131  -  151 

e 

P1IN,  P2IN  -  Input  parameters  ACCIN  53 

PITCH 

-  Calculated  in  V-GAMMA  160 

V 

PSI 

-  Calculated  BURN  116 

q 

GEOP 

-  Calculated  SELECT  27 

q1 

GEOA 

-  Calculated  SELECT  28 

r  (0) 

RREF 

-  Calculate  V-GAMMA  132 

R 

-  Initial  interim  value  dur¬ 
ing  pulse  burn  calculation 

RIP 

-  Initial  value  after  burn 
phase  of  burn  iteration 

R2P 

-  Final  value  after  coast 
phase  of  burn  iteration 

R_ 

RCUV(I) 

-  Calculated  in  RENTRY  37-41 

c 

Script  T 

T012 

-  Calculate  V-GAMMA  109 

At 

DTIN 

-  Input  parameter  ACCIN  47 

At 

DTN 

-  Calculated  in  SELECT  54 

n 

Atpc 

Same  as 

Atn 
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At. 


At(0,-e2) 


0 

0 

0 


C 

1 


1 


02  1 


0 


IN 


l 


DTL  -  Calculate  in  V-GAMMA  112 

Calculate  in  subroutine  TDRIFT 
Subroutine  TDRIFT  (0j,02,At) 

Same  as  M  pg  26 

THETAC  -  Calculated  in  RENTRY  31 

TIP  -  Initial  vehicle  location 

after  burn  phase  of  burn 
iteration 

Calculated  in  BURN  86 

T2P  -  Final  vehicle  location 

after  coast  phase  of  burn 
iteration 

Calculated  in  BURN  99 

Same  as  0  02*  only  for  a  time  segmented 

AV  calculation 


A0  {0j  At} 
V(0) 


AV 

AV, 

n 

AV 


N 


Wpc 
W!  Wa 

WIN,  WzN 
Y 


Calculated  in  subroutine  LDRIFT 
Subroutine  LDRIFT  (02,  AT,  Y) 

VREF  -  Calculated  V-GAMMA  133 

VP  -  Calculated  in  BURN  64 

V2P  -  Calculated  in  BURN  106 

VC  -  Calculated  in  RENTRY  33 

DVIN(I)  -  Input  parameter  ACCIN  47 

DVN  -  Calculate  at  TEMP  7  in  BURN 


DVSV(I) 

WFRE 

WjIN,  W2IN 
OWN 

Y1 IN,  Y2IN 


Calculate  in  BURN  58 

Calculated  in  RENTRY  73 

Input  parameters  ACCIN  42 

Calculated  in  SELECT  56 

Input  parameters  ACCIN  54 
Yaw 

Calculated  in  V-GAMMA  170 


COORDINATE  SYSTEMS 


n  n  n » 
i.  j,  k  f 


RREFUUO )  \ 

RREFUU(2)  >  Defined  in  V-GAMMA  131-140 
RREFUUO)  / 


{ 


El  SV 
E2SV 
E3SV 


Define  in  V-GAMMA  143-147 
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L6JSM*T200. 

ACCOUNT (30HNS0N* 200703 lS-JF*«OfVS* 9791) 
ATTACH (0L0PL*VGAMMA*IDaNSCHLGJ*C¥»4) 
REDUEST (NEWPt*  *PF  > 

UAOATE(P*N,L*0) 

FTNtA«f ,t=0r 
FlLtlTAPEl «8T=C*RT*S) 

LOStTTPILES»TAPEl) 

L«OtPU*77777r 
RETURN (NEWPL) 

RETURN (OLDPL) 

ATTACH (OLDPL ♦ SKETCH* I0=NSCNLGJ*C¥*2) 
UPDATE (P*L*0) 

REM1NQ ( TAPEl* 

RETURN (LGO) 

FTN  < A ♦ I *L®0 ) 

FILfc<TAPEl*BTsC*RT*S> 

L.OStTTPlLES»TAPEl) 

L«0i 

ooooooooooooeooooooooo 


\ 


Figure  Bl.  System  Dayflle  for  the  V-GAMMA 
and  SKETCH  Sample  Problems 
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